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Background 

Quantitative myocardial mapping allows assessment 
of focal and diffuse fibrosis in the myocardium, by sam- 
pling the Ti relaxation curve using inversion [1] or 
saturation recovery (SR) preparation [2] or a combination 
of both [3], followed by the acquisition of multiple 
images with different contrasts, which are subsequently 
fitted to a parametric equation pixel-wise to yield the 
maps. In myocardial Ti mapping, there is a degree of 
freedom in selecting which points on the relaxation curve 
are sampled. However, this topic has not been studied. In 
this study, we sought to develop an estimation theoretic 



framework for optimal selection of sampling points and 
characterized the variance of the corresponding esti- 
mator for sampling of the SR curve. 

Methods 

Based on the signal model, y^ = a (1-b exp(-Xk/Ti))+nk, 
and the least squares model, we derived the Fisher infor- 
mation matrix [4] . This was used to derive the Bayesian 
Cramer-Rao bound [4] for the variance of the Ti estima- 
tor for Ti values of interest between 950 and 1250 ms 
(~pre-contrast myocardium). The bound was evaluated 
for the SASHA sequence [2] which allows sampling 





uniformly distributed points 
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uniform 


theory std 
wrt. uniform 
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1457 ± 7.7 


69.5 


1456 ± 7.4 


48.4 


0.69 


0.71 
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1144± 14.5 


56.1 


1130±7.1 


41.1 


0.73 


0.76 
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1151 ±11.5 


53.3 


1155 ±8.6 


43.2 


0.81 


0.76 


4 


729 ± 10.3 


31.3 


724 ±2.1 


26.3 


0.84 


0.86 
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980 ± 11.2 
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981 ± 10.4 


25.2 


0.73 


0.78 
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823 ± 13.3 


29.9 


822 ±7.7 


24.2 


0.81 


0.83 


7 


1148 ±18.4 


53.0 


1144 ±8.3 


37.8 


0.71 


0.76 
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1130 ± 10.6 


56.1 


1137 ± 10.5 


45.2 


0.81 


0.76 
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963 ± 13.8 


50.0 


962 ± 6.2 


35.5 


0.71 


0.79 



Figure 1 Results of the phantom imaging over vials with values > 700 ms using the proposed and uniform sampling strategies, 
where each acquisition was repeated 5 times. The ratio of the standard deviation of the Ti estimator for each proposed sampling strategy 
and that of the uniform sampling strategy is reported as "standard deviation (std) with respect to (wrt) uniform." There is a gain in using the 
proposed point selection strategy, which is significantly different than 1 (P < 0.001). The values match those predicted by theory (P = 0.23). 
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uniform point selection 


proposed point selection 




subject 


anatomy 


— est , V 

Ti (ms) 


std(Ti''*) (ms) 


— est , . 

Ti (ms) 


std(Ti''*) (ms) 


std wrt. 
uniform 
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myocardium 


1211 ±6.4 


110.1 ± 10.3 


1201 ±9.2 


84.0 ± 8.2 


0 76 




blood 


1925 ±22.5 


166.7 ± 14.8 


1903 ± 17.2 


111.8±6.9 


0.67 
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myocardium 


1242 ±56.9 


132.9 ±32.5 


1254 ±58.9 


101.5 ± 12.4 


0.76 




blood 


1772 ±21.6 


204.9 ± 28.9 
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myocardium 


1187 ±55.0 
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1809 ±31.7 


138.1 ± 13.4 
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1168 ±48.5 
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76.6 ± 4.4 
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1772 ±42.7 


164.7 ± 16.9 


1761 ±20.8 


114.2 ±9.3 


0.69 



Figure 2 Results of in-vivo imaging for five healthy subjects using the proposed and uniform sampling strategies, where each 
acquisition was repeated 5 times. Ti^^^ is reported as the mean ± std of the average Ti values in the ROI across 5 scans, as a surrogate for 
accuracy and inter-scan reproducibility. The stddi^^^) is reported as the mean ± std of the std of the Ti values in the ROI across 5 scans, as a 
surrogate for the precision within the scan. Std wrt. uniform is the ratio of the mean values of stddi^^^) using the proposed and uniform point 
selection, as a surrogate for the percentage gain in precision. The standard deviation of the Ti estimator in the myocardium and blood was 
reduced by 23.6% and 26.8% respectively using the proposed approach. 



within a heart-beat between and T^ax with one 
point at full magnetization recovery (x^ = '^), and mini- 
mized over the choice of sampling points {xiJ yielding 
the proposed point selection. Phantom imaging of NiCl2 
doped agarose vials was performed to compare the pro- 
posed point selection with a uniform distribution of sam- 
pling points between T^^^ and T^ax [3] using an SSFP 
sequence with body-coil (NSA = 5) for 11 sampling 
points. Standard deviation (std) of Ti values within the 
vials was used as a surrogate for the variance of the esti- 
mator. Imaging was also performed on 5 healthy adult 
subjects (4 women, 23.4 ± 3.3 years) with a 32-channel 
cardiac-coil to verify the gains predicted by the theory. 
Both proposed and uniform point selection acquisitions 
were repeated 5 times per subject to average out the 
effects of noise. ROIs were drawn in the myocardium and 
the blood. Both the Ti estimate (average Ti values in the 
ROI) and the std of the estimator (std of values in the 
ROI) are reported as mean ± std across 5 scans. 

Results 

The point selection yielded a tri-modal distribution of 
points: 4 at T^^^, 6 at T^axj 1 at oo, with a theoretical 
gain in std of 24% compared to uniform selection. 
Figure 1 shows the results of phantom imaging for Ti 
values > 700 ms, indicating a good match between the- 
ory and experiment. Figure 2 depicts the measurements 
from the in-vivo data, averaged over five scans. Overall, 
there was a 23.6% and 26.8% reduction in the std of the 
Ti maps in the myocardium and blood respectively 
using the proposed approach. 



Conclusions 

The proposed framework allows for choosing the loca- 
tion of points on the Ti relaxation curve to achieve 
higher levels of precision without increasing the scan 
time. 
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